Experimental evaluation of the plunger technique: A method of cyclic manual aspiration thrombectomy for treatment of acute ischemic stroke.

BACKGROUND: Mechanical thrombectomy via direct aspiration is a rapid treatment for acute ischemic stroke. This method often results in the partial ingestion of the clot or "corking" of the catheter tip. Cyclic aspiration may take advantage of the mechanical properties of the clot, resulting in greater clot ingestion and overall procedure success. METHODS: An in vitro analysis was performed comparing static and cyclic (plunger technique) aspiration. Embolus analogs were used to create occlusions in a mock circulatory flow loop, and one aspiration attempt (first pass effect) using either a static or plunger technique was performed. The percent ingestion of each embolus analog was recorded for each trial. RESULTS: Static aspiration for 0% and 50% hematocrit embolus analogs resulted in ingestions of 12.8 ± 4.6% and 15.1 ± 10.0%, respectively, while plunger technique (cyclic) aspiration resulted in 15.8 ± 7.3% and 34.4 ± 19.5% ingestion. Complete ingestion was observed only with 50% hematocrit analogs, occurring in 30% of plunger and 10% of static cases. Statistical differences were determined between the two aspiration techniques for the 50% hematocrit samples, with the plunger technique yielding significantly more ingestion. In addition, the plunger technique was shown to maintain a negative vacuum pressure throughout the duration of cyclic plunging. CONCLUSIONS: The plunger technique for manual cyclic aspiration resulted in higher rates of complete ingestion and greater average % ingestions when compared to static aspiration. Increased clot ingestion can result in a higher rate of complete reperfusion during the first aspiration attempt, maximizing the number of patients with good clinical outcomes.

Analysis of frictional forces in experimental models of stent retriever mechanical thrombectomy.

Acute ischemic stroke (AIS) and mechanical thrombectomy (MT) are commonly studied in vitro using cerebral artery models made of nonbiological materials. However, these models often report higher recanalization rates than those observed clinically, suggesting a discrepancy between experimental models and clinical settings. We believe this may be partly due to the frictional interactions between blood clots, stent retrievers (SRs), and the vessel walls. Experiments were performed to measure the coefficients of static and kinetic friction between blood clots, common nonbiological model materials, and bovine carotid arteries (CAs). Additional friction testing was performed with nitinol SRs. Coefficients of static friction between blood clots and nonbiological materials were found to range from 0.1 to 0.2, increasing with decreasing clot hematocrit, but were significantly higher between blood clots and CAs (1.49, 0.57, and 0.46 for 0, 20, and 40% hematocrit clots, respectively). For 0% and 40% hematocrit clots, the coefficients of kinetic friction with nonbiological materials were less than 0.1, while significantly higher with CAs (0.26 and 0.23 for 0% and 40% hematocrit clots, respectively). However, no significant differences in the coefficients of kinetic friction were found between the different hematocrit clots. Testing with the nitinol SR showed no significant differences in the coefficients of kinetic friction for CAs (0.73) and silicone (0.78), suggesting that silicone could be a suitable model material for evaluating SR-vessel interactions in vitro. Overall, it is evident that discrepancies exist in the frictional forces between materials commonly used in experimental models of AIS and MT and those seen in vivo. The individual contributions of clot-artery, SR-artery, and clot-SR interactions during blood clot removal merit further investigation.

Effect of Hematocrit and Elevated Beat Rate on the 12cc Penn State Pediatric Ventricular Assist Device.

Congenital heart disease affects approximately 40,000 infants annually in the United States with 25% requiring invasive treatment. Due to limited number of donor hearts and treatment options available for children, pediatric ventricular assist devices (PVADs) are used as a bridge to transplant. The 12cc pneumatic Penn State PVAD is optimized to prevent platelet adhesion and thrombus formation at patient nominal conditions; however, children demonstrate variable blood hematocrit and elevated heart rates. Therefore, with pediatric patients exhibiting greater variability, particle image velocimetry is used to evaluate the PVAD with three non-Newtonian hematocrit blood analogs (20%, 40%, and 60%) and at two beat rates (75 and 120 bpm) to understand the device's performance. The flow fields demonstrate a strong inlet jet that transitions to a solid body rotation during diastole. During systole, the rotation dissipates and reorganizes into an outlet jet. This flow field is consistent across all hematocrits and beat rates but at a higher velocity magnitude during 120 bpm. There are also minor differences in flow field timing and surface washing due to hematocrit. Therefore, despite patient differences in hematocrit or required pumping output, thorough surface washing can be achieved in the PVAD by altering operating conditions, thus reducing platelet adhesion potential.

A MALDI Mass Spectrometry Imaging Sample Preparation Method for Venous Thrombosis with Initial Lipid Characterization of Lab-Made and Murine Clots.

Venous thromboembolism (VTE) and its complications affect over 900,000 people in the U.S. annually, with a third of cases resulting in fatality. Despite such a high incidence rate, venous thrombosis research has not led to significant changes in clinical treatments, with standard anti-coagulant therapy (heparin followed by a vitamin K antagonist) being used since the 1950s. Mechanical thrombectomy is an alternative strategy for treating venous thrombosis; however, clinical guidelines for patient selection have not been well-established or accepted. The effectiveness of both treatments is impacted by the heterogeneity of the thrombus, including the mechanical properties of its cellular components and its molecular makeup. A full understanding of the complex interplay between disease initiation and progression, biochemical molecular changes, tissue function, and mechanical properties calls for a multiplex and multiscale approach. In this work, we establish a protocol for using matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging to characterize spatial heterogeneity of biomolecules in lab-made blood clots and ex vivo murine thrombi. In this work, we compared (1) tissue preservation and cryosectioning methods, (2) various matrixes, 9-aminoacridine hydrochloride monohydrate (9AA), 2,5-dihydroxybenzoic acid (DHB), and alpha-cyano-4-hydroxycinnamic acid matrix (CHCA), (3) plasma-rich versus red-blood-cell rich lab-made blood clots, and (4) lab-made blood clots versus ex vivo murine thrombi. This project is the first step in our work to combine mass spectrometry imaging with biomechanical testing of blood clots to improve our understanding of VTE.

The influence of blood composition and loading frequency on the behavior of embolus analogs.

In cases of acute ischemic stroke (AIS), mechanical thrombectomy (MT) can be used to directly remove lodged thromboemboli. Despite improvements in patient outcomes, one of the key factors affecting MT success is the mechanical properties of the occlusive thrombus. Therefore, the goal of this study was to investigate the viscoelastic properties of embolus analogs (EAs) and determine the influence of EA hematocrit and loading frequency. Bovine blood EAs were created over a range of physiological hematocrits (0-60%) and cyclic uniaxial compression testing was performed at three loading frequencies to mimic in vivo loading conditions, followed by stress-relaxation testing. It was found that EAs exhibited behaviors typical of hyper-viscoelastic materials and that EA hematocrit played a large role in both EA stiffness and relaxation, with both parameters decreasing as hematocrit increased from 0 to 60%. The viscoelastic behavior of the EAs was also affected by the frequency at which they were loaded, with significant increases in peak stresses between the 0.5 and 2 Hz loaded EAs. Lower hematocrit EAs had very dense fibrin networks while the higher hematocrit EAs consisted of closely packed RBCs with little fibrin present. These results suggest that fibrin contributes to EA stiffness and relaxation behaviors while RBCs play a role in decreasing the overall viscous response and strain-rate dependency. An Ogden hyperelastic model was found to best reproduce the EA loading data while a 3-term Prony series was fit to the stress relaxation data. A hyper-viscoelastic modeling framework was then implemented combining the loading and stress-relaxation fits and the results could match the full cyclic loading data for EAs of varying hematocrit and loading frequency. The results of the experimental mechanical characterization and hyper-viscoelastic curve fitting can be incorporated in future modeling efforts to optimize mechanical thrombectomy for AIS patients.

The effects of non-Newtonian blood modeling and pulsatility on hemodynamics in the food and drug administration's benchmark nozzle model.

BACKGROUND: Computational fluid dynamics (CFD) is an important tool for predicting cardiovascular device performance. The FDA developed a benchmark nozzle model in which experimental and CFD data were compared, however, the studies were limited by steady flows and Newtonian models. OBJECTIVE: Newtonian and non-Newtonian blood models will be compared under steady and pulsatile flows to evaluate their influence on hemodynamics in the FDA nozzle. METHODS: CFD simulations were validated against the FDA data for steady flow with a Newtonian model. Further simulations were performed using Newtonian and non-Newtonian models under both steady and pulsatile flows. RESULTS: CFD results were within the experimental standard deviations at nearly all locations and Reynolds numbers. The model differences were most evident at Re = 500, in the recirculation regions, and during diastole. The non-Newtonian model predicted blunter upstream velocity profiles, higher velocities in the throat, and differences in the recirculation flow patterns. The non-Newtonian model also predicted a greater pressure drop at Re = 500 with minimal differences observed at higher Reynolds numbers. CONCLUSIONS: An improved modeling framework and validation procedure were used to further investigate hemodynamics in geometries relevant to cardiovascular devices and found that accounting for blood's non-Newtonian and pulsatile behavior can lead to large differences in predictions in hemodynamic parameters.

Computational Investigation of Anastomosis Options of a Right-Heart Pump to Patient Specific Pulmonary Arteries.

Patients with Fontan circulation have increased risk of heart failure, but are not always candidates for heart transplant, leading to the development of the subpulmonic Penn State Fontan Circulation Assist Device. The aim of this study was to use patient-specific computational fluid dynamics simulations to evaluate anastomosis options for implanting this device. Simulations were performed of the pre-surgical anatomy as well as four surgical options: a T-junction and three Y-grafts. Cases were evaluated based on several fluid-dynamic quantities. The impact of imbalanced left-right pulmonary flow distribution was also investigated. Results showed that a 12-mm Y-graft was the most energy efficient. However, an 8-mm graft showed more favorable wall shear stress distribution, indicating lower risk of thrombosis and endothelial damage. The 8-mm Y-grafts also showed a more balanced pulmonary flow split, and lower residence time, also indicating lower thrombosis risk. The relative performance of the surgical options was largely unchanged whether or not the pulmonary vascular resistance remained imbalanced post-implantation.

Computational Modeling of the Penn State Fontan Circulation Assist Device.

To address the increasing number of failing Fontan patients, Penn State University and the Penn State Hershey Medical Center are developing a centrifugal blood pump for long-term mechanical support. Computational fluid dynamics (CFD) modeling of the Penn State Fontan Circulatory Assist Device (FCAD) was performed to understand hemodynamics within the pump and its potential for hemolysis and thrombosis. CFD velocity and pressure results were first validated against experimental data and found to be within the standard deviations of the velocities and within 5% of the pressures. Further simulations performed with a human blood model found that most of the fluid domain was subjected to low shear stress (<50 Pa), with areas of highest stress around the rotor blade tips that increased with pump flow rate and rotor speed (138-178 Pa). However, the stresses compared well to previous CFD studies of commercial blood pumps and remained mostly below common thresholds of hemolysis and platelet activation. Additionally, few regions of low shear rate were observed within the FCAD, signifying minimal potential for platelet adhesion. These results further emphasize the FCAD's potential that has been observed previously in experimental and animal studies.

Experimental Hemodynamics Within the Penn State Fontan Circulatory Assist Device.

For children born with a single functional ventricle, the Fontan operation bypasses the right ventricle by forming a four-way total cavopulmonary connection and adapts the existing ventricle for the systemic circulation. However, upon reaching adulthood, many Fontan patients exhibit low cardiac output and elevated venous pressure, eventually requiring a heart transplantation. Despite efforts in developing a new device or using an existing device for failing Fontan support, there is still no Food and Drug Administration-approved device for subpulmonary support. Penn State University is developing a hydrodynamically levitated Fontan circulatory assist device (FCAD) for bridge-to-transplant or destination therapy. The hemodynamics within the FCAD, at both steady and patient averaged pulsatile conditions for three physiological pump operating conditions, were quantified using particle image velocimetry (PIV) to determine the velocity magnitudes and Reynolds normal and shear stresses within the device. Data were acquired at three planes (0 mm and ±25% of the radius) for the inferior and superior vena cavae inlets and the pulmonary artery outlet. The inlets had a blunt velocity profile that became skewed toward the collecting volute as fluid approached the rotor. At the outlet, regardless of the flow condition, a high-velocity jet exited the volute and moved downstream in a helical pattern. Turbulent stresses observed at the volute exit were influenced by the rotor's rotation. Regardless of inlet conditions, the pump demonstrated advantageous behavior for clinical use with a predictable flow field and a low risk of platelet adhesion and hemolysis based on calculated wall shear rates and turbulent stresses, respectively.

Hydrodynamics in Acute Ischemic Stroke Catheters Under Static and Cyclic Aspiration Conditions.

PURPOSE: Previous studies have suggested improved recanalization efficiency by using cyclic aspiration as opposed to static aspiration for mechanical thrombectomy in the treatment of acute ischemic stroke (AIS). However, there has not been an investigation into which parameters optimize this theoretical approach. Therefore, this study was designed to investigate the application of static and cyclic pressure in AIS aspiration catheters. METHODS: An experimental flow system was designed to apply and measure both static and cyclic aspiration pressures in four commercial aspiration catheters. A medical aspirator and solenoid-valve were used to generate a range of applied pump pressures and cyclic valve frequencies. The measured mean and oscillatory pressure amplitudes were analyzed and compared across catheters and operating conditions. RESULTS: The greatest mean pressures and forces were achieved under contact-aspiration conditions using the largest diameter catheter, while the largest oscillating pressure (ΔP) and force (ΔF) values were generated in the least compliant catheter, determined to be proportional to both its inner radius and wall thickness. Furthermore, lower frequency cyclic aspiration was able to generate the largest values of ΔP and ΔF. CONCLUSION: For future investigations of cyclic aspiration AIS therapy, the best options are to use large inner diameter and low compliance catheters operating at lower cyclic frequencies to generate the greatest mean and oscillating forces on lodged thromboemboli.

Computational modeling of the Food and Drug Administration's benchmark centrifugal blood pump.

In order to simulate hemodynamics within centrifugal blood pumps and to predict pump hemolysis, CFD simulations must be thoroughly validated against experimental data. They must also account for and accurately model the specific working fluid in the pump, whether that is a blood-analog solution to match an experimental PIV study or animal blood in a hemolysis experiment. Therefore, the Food and Drug Administration (FDA) benchmark centrifugal blood pump and its database of experimental PIV and hemolysis data were used to thoroughly validate CFD simulations of the same blood pump. A Newtonian blood model was first used to compare to the PIV data with a blood analog fluid while hemolysis data were compared using a power-law hemolysis model fit to porcine blood data. A viscoelastic blood model was then incorporated into the CFD solver to investigate the importance of modeling blood's viscoelasticity in centrifugal pumps. The established computational framework, including a dynamic rotating mesh, animal blood-specific fluid properties and hemolysis modeling, and a k-ω SST turbulence model, was shown to more accurately predict pump pressure heads, velocity fields, and hemolysis compared to previously published CFD studies of the FDA centrifugal pump. The CFD simulations were able to match the FDA pressure and hemolysis data for multiple pump operating conditions, with the CFD results being within the standard deviations of the experimental results. While CFD radial velocity profiles between the impeller blades also compared well to the PIV velocity results, more work is still needed to address the large variability among both experimental and computational predictions of velocity in the diffuser outlet jet. Small differences were observed between the Newtonian and viscoelastic blood models in pressure head and hemolysis at the higher flow rate cases (FDA Conditions 4 and 5) but were more significant at lower flow rate and pump impeller speeds (FDA Condition 1). These results suggest that the importance of accounting for blood's viscoelasticity may be dependent on the specific blood pump operating conditions. This detailed computational framework with improved modeling techniques and an extensive validation procedure will be used in future CFD studies of centrifugal blood pumps to aid in device design and predictions of their biological responses.

Development of a computational model for acute ischemic stroke recanalization through cyclic aspiration.

Acute ischemic stroke (AIS), the result of embolic occlusion of a cerebral artery, is responsible for 87% of the 6.5 million stroke-related deaths each year. Despite improvements from first-generation thrombectomy devices for treating AIS, 80% of eligible stroke patients will either die or suffer a major disability. In order to maximize the number of patients with good outcomes, new AIS therapies need to be developed to achieve complete reperfusion on the first pass. One such therapy that has shown promise experimentally is the application of cyclic aspiration pressure, which led to higher recanalization rates at lower pressure magnitudes. In order to investigate AIS and cyclic aspiration recanalization, an improved computational modeling framework was developed, combining a viscoelastic thromboembolus model with a cohesive zone (CZ) model for the thromboembolus-artery interface. The model was first validated against experimental displacement data of a cyclically aspirated thromboembolus analog. The CZ model parameters, including the addition of a damage accumulation model, were then investigated computationally to determine their individual effects on the thromboembolus and CZ behavior. The relaxation time and the damage model critical opening length were shown to have the greatest effect on the CZ opening and led to increased displacement that accumulated with repeated loading. Additional simulations were performed with parameters relevant to AIS including internal carotid artery dimensions and thromboemboli mechanical properties. In these AIS cases, more upstream CZ opening was observed compared to the thromboembolus analog cases and greater displacement was achieved with the lower-frequency aspiration (0.5 vs 1 Hz).

Inter-Laboratory Characterization of the Velocity Field in the FDA Blood Pump Model Using Particle Image Velocimetry (PIV).

PURPOSE: A credible computational fluid dynamics (CFD) model can play a meaningful role in evaluating the safety and performance of medical devices. A key step towards establishing model credibility is to first validate CFD models with benchmark experimental datasets to minimize model-form errors before applying the credibility assessment process to more complex medical devices. However, validation studies to establish benchmark datasets can be cost prohibitive and difficult to perform. The goal of this initiative sponsored by the U.S. Food and Drug Administration is to generate validation data for a simplified centrifugal pump that mimics blood flow characteristics commonly observed in ventricular assist devices. METHODS: The centrifugal blood pump model was made from clear acrylic and included an impeller, with four equally spaced, straight blades, supported by mechanical bearings. Particle Image Velocimetry (PIV) measurements were performed at several locations throughout the pump by three independent laboratories. A standard protocol was developed for the experiments to ensure that the flow conditions were comparable and to minimize systematic errors during PIV image acquisition and processing. Velocity fields were extracted at the pump entrance, blade passage area, back gap region, and at the outlet diffuser regions. A Newtonian blood analog fluid composed of sodium iodide, glycerin, and water was used as the working fluid. Velocity measurements were made for six different pump flow conditions, with the blood-equivalent flow rate ranging between 2.5 and 7 L/min for pump speeds of 2500 and 3500 rpm. RESULTS: Mean intra- and inter-laboratory variabilities in velocity were ~ 10% at the majority of the measurement locations inside the pump. However, the inter-laboratory variability increased to more than ~ 30% in the exit diffuser region. The variability between the three laboratories for the peak velocity magnitude in the diffuser region ranged from 5 to 25%. The bulk velocity field near the impeller changed proportionally with the rotational speed but was relatively unaffected by the pump flow rate. In contrast, flow in the exit diffuser region was sensitive to both the flow rate and the rotational speed. Specifically, at 3500 rpm, the exit jet tilted toward the inner wall of the diffuser at a flow rate of 2.5 L/min, but the jet tilted towards the outer wall when the flow rate was 7 L/min. CONCLUSIONS: Inter-laboratory experimental mean velocity data (and the corresponding variance) were obtained for the FDA pump model and are available for download at https://nciphub.org/wiki/FDA_CFD . Experimental datasets from the inter-laboratory characterization of benchmark flow models, including the blood pump model presented herein and our previous nozzle model, can be used for validating future CFD studies and to collaboratively develop guidelines on best practices for verification, validation, uncertainty quantification, and credibility assessment of CFD simulations in the evaluation of medical devices (e.g. ASME V&V 40 standards working group).

Characterizing the HeartMate II Left Ventricular Assist Device Outflow Using Particle Image Velocimetry.

Ventricular assist devices (VADs) are implanted in patients with a diseased ventricle to maintain peripheral perfusion as a bridge-to-transplant or as destination therapy. However, some patients with continuous flow VADs (e.g., HeartMate II (HMII)) have experienced gastrointestinal (GI) bleeding, in part caused by the proteolytic cleavage or mechanical destruction of von Willebrand factor (vWF), a clotting glycoprotein. in vitro studies were performed to measure the flow located within the HMII outlet cannula under both steady and physiological conditions using particle image velocimetry (PIV). Under steady flow, a mock flow loop was used with the HMII producing a flow rate of 3.2 L/min. The physiological experiment included a pulsatile pump operated at 105 BPM with a ventricle filling volume of 50 mL and in conjunction with the HMII producing a total flow rate of 5.0 L/min. Velocity fields, Reynolds normal stresses (RNSs), and Reynolds shear stresses (RSSs) were analyzed to quantify the outlet flow's potential contribution to vWF degradation. Under both flow conditions, the HMII generated principal Reynolds stresses that are, at times, orders of magnitude higher than those needed to unfurl vWF, potentially impacting its physiological function. Under steady flow, principal RNSs were calculated to be approximately 500 Pa in the outlet cannula. Elevated Reynolds stresses were observed throughout every phase of the cardiac cycle under physiological flow with principal RNSs approaching 1500 Pa during peak systole. Prolonged exposure to these conditions may lead to acquired von Willebrand syndrome (AvWS), which is accompanied by uncontrollable bleeding episodes.

Asynchronous Pumping of a Pulsatile Ventricular Assist Device in a Pediatric Anastomosis Model.

BACKGROUND: Both pulsatile and continuous flow ventricular assist devices are being developed for pediatric congenital heart defect patients. Pulsatile devices are often operated asynchronously with the heart in either an "automatic" or a fixed beat rate mode. However, most studies have only investigated synchronized ejection. METHODS: A previously validated viscoelastic blood solver is used to investigate the parameters of pulsatility, power loss, and graft failure in a pediatric aortic anastomosis model. RESULTS: Pulsatility was highest with synchronized flow and lowest at a 90° phase shift. Power loss decreased at 90° and 180° phase shifts but increased at a 270° phase shift. Similar regions of potential intimal hyperplasia and graft failure were seen in all cases but with phase-shifted ejection leading to higher wall shear stress on the anastomotic floor and oscillatory shear index on the anastomotic toe. CONCLUSION: The ranges of pulsatility and hemodynamics that can result clinically using asynchronous pulsatile devices were investigated in a pediatric anastomosis model. These results, along with the different postoperative benefits of pump modulation, can be used to design an optimal weaning protocol.

FDA Benchmark Medical Device Flow Models for CFD Validation.

Computational fluid dynamics (CFD) is increasingly being used to develop blood-contacting medical devices. However, the lack of standardized methods for validating CFD simulations and blood damage predictions limits its use in the safety evaluation of devices. Through a U.S. Food and Drug Administration (FDA) initiative, two benchmark models of typical device flow geometries (nozzle and centrifugal blood pump) were tested in multiple laboratories to provide experimental velocities, pressures, and hemolysis data to support CFD validation. In addition, computational simulations were performed by more than 20 independent groups to assess current CFD techniques. The primary goal of this article is to summarize the FDA initiative and to report recent findings from the benchmark blood pump model study. Discrepancies between CFD predicted velocities and those measured using particle image velocimetry most often occurred in regions of flow separation (e.g., downstream of the nozzle throat, and in the pump exit diffuser). For the six pump test conditions, 57% of the CFD predictions of pressure head were within one standard deviation of the mean measured values. Notably, only 37% of all CFD submissions contained hemolysis predictions. This project aided in the development of an FDA Guidance Document on factors to consider when reporting computational studies in medical device regulatory submissions. There is an accompanying podcast available for this article. Please visit the journal's Web site (www.asaiojournal.com) to listen.

Analysis of Transitional and Turbulent Flow Through the FDA Benchmark Nozzle Model Using Laser Doppler Velocimetry.

Transitional and turbulent flow through a simplified medical device model is analyzed as part of the FDA's Critical Path Initiative, designed to improve the process of bringing medical products to market. Computational predictions are often used in the development of devices and reliable in vitro data is needed to validate computational results, particularly estimations of the Reynolds stresses that could play a role in damaging blood elements. The high spatial resolution of laser Doppler velocimetry (LDV) is used to collect two component velocity data within the FDA benchmark nozzle model. Two flow conditions are used to produce flow encompassing laminar, transitional, and turbulent regimes, and viscous stresses, principal Reynolds stresses, and turbulence intensities are calculated from the measured LDV velocities. Axial velocities and viscous stresses are compared to data from a prior inter-laboratory study conducted with particle image velocimetry. Large velocity gradients are observed near the wall in the nozzle throat and in the jet shear layer located in the expansion downstream of the throat, with axial velocity changing as much as 4.5 m/s over 200 µm. Additionally, maximum Reynolds shear stresses of 1000-2000 Pa are calculated in the high shear regions, which are an order of magnitude higher than the peak viscous shear stresses (<100 Pa). It is important to consider the effects of both viscous and turbulent stresses when simulating flow through medical devices. Reynolds stresses above commonly accepted hemolysis thresholds are measured in the nozzle model, indicating that hemolysis may occur under certain flow conditions. As such, the presented turbulence quantities from LDV, which are also available for download at https://fdacfd.nci.nih.gov/ , provide an ideal validation test for computational simulations that seek to characterize the flow field and to predict hemolysis within the FDA nozzle geometry.

Hemodynamics in a Pediatric Ascending Aorta Using a Viscoelastic Pediatric Blood Model.

Congenital heart disease is the leading cause of infant death in the United States with over 36,000 newborns affected each year. Despite this growing problem there are few mechanical circulatory support devices designed specifically for pediatric and neonate patients. Previous research has been done investigating pediatric ventricular assist devices (PVADs) assuming blood to be a Newtonian fluid in computational fluid dynamics (CFD) simulations, ignoring its viscoelastic and shear-thinning properties. In contrast to adult VADs, PVADs may be more susceptible to hemolysis and thrombosis due to altered flow into the aorta, and therefore, a more accurate blood model should be used. A CFD solver that incorporates a modified Oldroyd-B model designed specifically for pediatric blood is used to investigate important hemodynamic parameters in a pediatric aortic model under pulsatile flow conditions. These results are compared to Newtonian blood simulations at three physiological pediatric hematocrits. Minor differences are seen in both velocity and wall shear stress (WSS) during early stages of the cardiac cycle between the Newtonian and viscoelastic models. During diastole, significant differences are seen in the velocities in the descending aorta (up to 12%) and in the aortic branches (up to 30%) between the two models. Additionally, peak WSS differences are seen between the models throughout the cardiac cycle. At the onset of diastole, peak WSS differences of 43% are seen between the Newtonian and viscoelastic model and between the 20 and 60% hematocrit viscoelastic models at peak systole of 41%.